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O ' Abstract 

o : 

■ Triple-gluon elastic scatterings are briefly reviewed since the scat- 

terings explain the early thermalization puzzle in Au-Au collisions at 
RHIC energies. A numerical solution of the transport equation with the 
triple-gluon elastic scatterings demonstrates gluon momentum isotropy 
achieved at a time of the order of 0.65 fm/c. Triple-gluon scatterings 
, lead to a short thermalization time of gluon matter. 
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Hadron spectra at low momentum can be fitted to thermal distributions. 
The disappearance of angular correlation of back-to-back jets means that one 
jet is dissolved into medium and thermalized. The thermalization of matter 
created in Au-Au collisions at RHIC is obvious. However, relevant importance 
is that the initially created matter thermalizes within a time of less than 1 
fm/c [1,2] at least in the regime near midrapidity [3]. The early thermaliza- 
tion is concluded while the results of hydrodynamic calculations [1,3-7] are 
in agreement with the elliptic flow data [8,9] at p± < 2 GeV/c. Triple-gluon 
elastic scatterings were proposed to explain the early thermalization puzzle in 
Ref. [10]. 

Why the gluon-gluon-gluon elastic scattering processes are needed at 
RHIC energies? The three-gluon to three-gluon scatterings get important while 
the gluon number density is high. To account for RHIC experimental data 
the gluon rapidity density dN s /dy ~ 1000 is suggested for the initial gluon 
matter [11,12]. Such a value of the rapidity density leads to a gluon number 
density of about 38 fm -3 , which is high enough for triple-gluon scatterings to 
occur. When the gluon number density gets larger, three-gluon scatterings 
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become more important. The three-gluon scattering processes are anticipated 
to thoroughly overwhelm the two-gluon elastic scattering processes in heavy ion 
collisions at LHC energies. An estimate from minijet production gives initial 
gluon number densities of about 30 fm -3 for RHIC and 140 fm~ 3 for LHC [13]. 
This supports the importance of the three-gluon scattering processes. 
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Fig. 1: Scatterings of three gluons. 

2. Triple-gluon scatterings 

The two-gluon to two-gluon scattering processes were studied in per- 
turbative QCD by Cutler and Sivers [14], Combridge, Kripfganz and Ranft 
[15]. The spin- and color-averaged squared amplitude | | 2 for the two- 

gluon elastic scatterings at order a 2 s was expressd in terms of the Mandel- 
stam variables. For the triple-gluon elastic scatterings, many diagrams have 



2 



to be calculated by fortran codes. Some of the three-gluon to three-gluon 
scattering diagrams are presented in Figs. 1 and 2 and corresponding spin- 
and color- averaged squared amplitude | ^3^3 | 2 is calculated in perturba- 
tive QCD. A gluon four-momentum is labeled as pi = (Ei,p{) in the process 
g(Pi) + g(P2) + gfe) -> g(P4) + g(ps) + g(Pe)- Then nine Lorentz-invariant 
variables are defined as 

S12 = (Pi + P2f, S 2 3 = (P2 + P3f, S31 = (p 3 + Pl) 2 
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These independent variables are appropriate for expressing | .M3-+3 
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Fig. 2: Scatterings of three gluons. 

Only the diagram in Fig. 2 is calculable by hand. This process 
consists of two four-gluon couplings and one gluon propagator between the 
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two vertices. The spin- and color-averaged squared amplitude | 1 2 for 

the diagram then depends on the three variables s i2 , s 23 and s 3 i in the 
simple form 

,J ' B ~~ 1 128 (s 12 + S23 + S31 ) 2 
which is independent of the momenta of the three final gluons. This ensures 
that each of the final gluons runs away in any direction in momentum space 
with equal probability. In other words, local momentum isotropy must be 
attained through the scattering process. Here, g s is the quark-gluon coupling 
constant, g 2 = 4na s and a s = 0.3 is used 

3. Thermalization of gluon matter 

Since the triple-gluon elastic scatterings are important in the initial gluon 
matter, a transport equation of Boltzmann type must include the 3-^3 scat- 
tering term, 
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2E lg33 J (2n) 3 2E 2 {2n) 3 2E 3 (2n) 3 2E 4 (2n) 3 2E 5 (2n) 3 2E 6 
x (2n) 4 S 4 (p 1 + p 2 + p 3 - p 4 - p 5 - p 6 ) | M 3 ^ 3 | 2 
x [/i/ 2 / s (l + / 4 )(1 + / 5 )(1 + /e) - /4/ 5 / 6 (l + /i)(l + / 2 )(1 + / 3 )] 

(2) 

where g 22 = 2, g 33 = 12, the degeneracy factor g G = 16 and the velocity of 
a massless gluon vi = 1. The gluon distribution function / ; depends on the 
position n, the momentum p\ and the time t. The first term on the right-hand- 
side of the equation exhibits the well-known 2^2 scattering processes. The 
second term is a new term, which represents the 3 — > 3 scatterings. Six gluon 
distributions are involved in the new term. The 3^3 scattering processes 
involve a larger phase space than the 2-^2 scattering processes. 

The squared amplitude | .M 3 ^ 3 | 2 has a very much long expression and 
substantially enhance the running time of computer in solving the transport 
equation. To gain a first sight at the physics determined by the transport 
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equation, initially produced gluons are assumed to uniformly distribute in 
a cylinder formed in a central heavy ion collision. Consequently, the gluon 
distribution function is only a function of momentum magnitude and time. We 
solve the equation starting at the time t ni \ when gluon matter is created and 
ending at the time U so when gluon matter reaches local momentum isotropy. 
The transverse expansion of the initially produced gluon matter is neglected 
due to a short thermalization time. 




Fig. 3: Gluon distribution functions versus momentum in different directions 
while gluon matter is just produced in the initial Au-Au collision. The dot- 
ted, dashed and dot-dashed curves correspond to the angles 6 = 0°, 45°, 90°, 
respectively. The solid curve stands for the thermal distribution function. 

For a central Au-Au collision at ^/snn = 200 GeV, an initial gluon 
distribution for the transport equation is obtained from Eq. (34) of Ref. [16], 

mtbd) = M * 10W- 5 e -b1/(0.9co sh ( y ))-(|p1tan h ( y ))V8g (y 2 _ 2) 

g G 7iR A Y(\ P | /cosh(y) + 0.3) 

(3) 

where | p | is in GeV, the nuclear radius Ra = 6.4 fm and y is the rapidity. 
The time t ini is 0.2 fm/c as estimated in HIJING Monte Carlo simulation [17]. 
Here, Y is the maximum of rapidity and approximately equals 5. 6{Y 2 — y 2 ) 
is the step function which is zero for Y < y or 1 for Y > y. Eq. (3) ex- 
hibits one form of initial gluon distribution in gluon matter which is not in 
thermal and chemical equilibrium. We draw in Fig. 3 the dotted, dashed and 
dot-dashed curves individually for the gluon distribution functions against the 
gluon momentum at the three angles 9 = 0°, 45°, 90° relative to the incom- 
ing beam direction. These curves can not coincide and do not exhibit local 
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momentum isotropy. Particularly, the gluon momentum distribution in the 
transverse direction differs considerably from one in the longitudinal direction. 
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Fig. 4: Gluon distribution functions versus momentum in different directions 
while gluon matter just arrives at thermal equilibrium. The dotted, dashed 
and dot-dashed curves correspond to the angles 6 = 0°,45 o ,90°, respectively. 
The solid curve represents the thermal distribution function. 

The difference among the initial gluon distribution functions in different 
directions does not last for a long time since the transport equation alters them. 
Momentum isotropy can be established at the time ti SO = 0.65 fm/c as shown 
by the dotted, dashed and dot-dashed curves in Fig. 4. The gluon distribution 
functions can be well fitted to the Juttner distribution with nonequilibrium 
fugacity A, 



where T gives the temperature of gluon matter. The solid curve in Fig. 4 
gives a fit of A = 0.065 and T = 0.75 GeV. We get the thermalization time 
^iso _ ^ini = 0.45 fm/c. The fugacity and the temperature are higher than 0.05 
and 0.55 GeV, and the thermalization time is less than 0.5 fm/c obtained from 
the free streaming of partons in Ref. [16], respectively. 

4. Summary 

The three-gluon scattering processes proposed to study the early ther- 
malization problem at high density have been reviewed. The gauge-invariant 
squared amplitude | M.^?, | 2 for all the triple-gluon scatterings is given at 
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the lowest order ag. | ^3^3 | 2 enters the transport equation to give a new 
contribution to the time dependence of the gluon distribution function. The 
evolution of gluon matter is dominated by the three-gluon scattering processes 
while the number density is high. With the initial gluon matter obtained from 
HIJING simulation, the transport equation gives the thermalization time of 
about 0.45 fm/c. The three-gluon scatterings considerably shorten the ther- 
malization time of gluon matter at high density. 
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